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OPTIMUM TRAJECTORIES TO CIRCULAR SYNCHRONOUS EQUATORIAL ORBIT 
FOR SMALLER-THAN-OPTIMUM APOGEE MOTORS 
byOmerF. Spurlock and Fred Teren 
Lewis Research Center 

SUMMARY 

Analysis, procedure, and results are presented for maximizing the payload capa- 
bility for trajectories to circular synchronous equatorial orbit where the apogee motor 
total impulse is much smaller than optimum. The trajectories begin at launch and are 
numerically integrated to insertion into the final orbit. Constraints on parking orbit 
perigee radius and duration are included as part of the solution. The trajectories, unlike 
conventional synchronous orbit trajectories, were found to require noncircular parking 
orbits and large amounts of inclination reduction before the apogee burn. A calculus of 
variations formulation of the optimization problem is used to obtain the solution. 

Results are presented for the Applications Technology Satellite (ATS)-E mission, a 
circular synchronous equatorial orbit mission. The total impulse of the solid apogee 
motor is much smaller than optimum. The spacecraft mass for this mission was fixed. 

A conventional trajectory could not place the spacecraft in the required orbit. The un- 
conventional trajectory was determined to be able to achieve the orbit. 

The analysis and procedure were also applied to the conventional circular synchro- 
nous equatorial orbit problem where the burn and coast durations are optimum. Results 
are presented for a mission of this type to show that the familiar conventional circular 
synchronous equatorial orbit launch profile is essentially optimum. 


INTRODUCTION 

A conventional trajectory to circular synchronous equatorial orbit consists of five 
consecutive phases as shown in figure 1. The first phase is an ascent from the launch 
site to a circular parking orbit. To maximize the mass in parking orbit, a 90° launch 
azimuth is used, which results in a parking orbit inclination equal to the launch site 
latitude. This inclination must be removed during the trajectory. The second phase 
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Figure 1. -Conventional trajectory to circular synchronous equatorial 
orbit. (Plane changes are not shown. ) 

is a coast arc to the proximity of the equator. A small portion of the required plane 
change is removed by the second burn, the third phase. Much more importantly, the 
second burn must place the vehicle in a transfer orbit whose apogee is over the equator 
and equal to synchronous altitude. The vehicle coasts to apogee of the transfer orbit in 
the fourth phase. The fifth phase consists of a final burn that removes the major portion 
of the inclination and circularizes the orbit. Intuitively, it seems reasonable that this 
conventional profile is near optimum if the burn and coast durations may be varied to 
maximize the mass at the end of each burn. However, if the total impulse of the final 
burn is fixed at less than the optimum value, the conventional trajectory must be mod- 
ified to yield maximum payload to the final orbit. In particular, the parking orbit is 
noncircular, the perigee radius of the transfer orbit increases, and the second burn 
removes more than a minor part of the inclination. The optimization problem is to find 
the best combination of these changes and other less important ones to yield maximum 
payload to circular synchronous equatorial orbit. 

Optimization of the conventional trajectory to circular synchronous equatorial orbit 
has been treated by several authors. Hoelker and Silber (ref. 1) present a detailed 
analysis of the conventional problem. Rider (ref. 2) considers the problem of changing 
the plane and also the radius of a circular orbit. These and other similar studies treat 
the problem as one of changing the plane and radius of a circular orbit, ignoring the 
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ascent to the first (parking) circular orbit. This is satisfactory for the conventional 
case. However, an unconventional trajectory is more complex since the parking orbit 
is in general noncircular. The ascent must be included as part of the optimization prob- 
lem. Therefore, a more sophisticated optimization procedure is required for uncon- 
ventional trajectories. Additionally, references 1 and 2 are general and consequently 
are not concerned with constraints which may alter the acceptability of a given trajec- 
tory, such as limitations on coast time or the minimum perigee radius of the noncircular 
parking orbit. 

The problem of optimizing trajectories to circular synchronous equatorial orbit may 
be considered as a multistage launch vehicle optimization in which two or more of the 
burns and coasts are free for optimization. Several analyses have been performed to 
optimize multistage launch vehicles, including one by the authors of this report (ref. 3). 
For optimizing the unconventional trajectory, the analysis in reference 3 was expanded 
to three dimensions and, also, to include a perigee radius constraint on the parking orbit. 
The perigee radius constraint must be included to limit aerodynamic heating on the ve- 
hicle. 

An example of a mission which requires an unconventional trajectory is the Appli- 
cations Technology Satellite (ATS)-E mission on the Atlas -Centaur vehicle. The final 
burn is performed by a solid motor which is part of the spacecraft system. That motor 
is significantly smaller than optimum. There are spacecraft and launch vehicle con- 
straints which must be incorporated into the solution. The perigee radius and the park- 
ing orbit coast duration are limited. The results for this problem are presented as an 
example of trajectory optimization for unconventional profiles. 

The more complex analysis was also applied to a case where the apogee motor is of 
optimum size. References 1 and 2 show that the conventional trajectory is optimum if the 
ascent to parking orbit is ignored and if all burns are treated impulsively. This report 
examines the case where these simplifications are not made. An optimum trajectory is 
obtained and described. 


ANALYSIS 

Problem Description 

A conventional trajectory to circular synchronous equatorial orbit consists of five 
phases: 

(1) Ascent to parking orbit 

(2) Parking orbit coast 

(3) Second impulse 
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Figure 2. - Circular synchronous equatorial orbit ascent profile. 


(4) Transfer orbit coast 

(5) Third impulse or apogee burn 

Figure 1 shows the planar characteristics of the conventional trajectory. The 
nonplanar characteristics are shown in figure 2. 

The vehicle is launched at an azimuth of 90° in order to maximize the vehicle mass 
in parking orbit and to minimize the inclination of the parking orbit. The circular park- 
ing orbit altitude is as low as aerodynamic heating constraints will allow, usually about 
165 to 185 kilometers. The parking orbit coast time is usually 15 minutes - the time 
required to coast from orbit insertion to the first equator crossing. The third phase 
places the vehicle in a transfer orbit whose apogee and perigee are over the equator. 

The apogee altitude is about equal to the required altitude for a circular synchronous 
orbit. The transfer orbit coast time is about 5^ hours. The third impulse, the apogee 
burn, occurs at apogee of the transfer orbit. Apogee is designed to occur at the equator 
and at the proper altitude for injection into the final orbit. A small part of the inclin- 
ation is removed by the second impulse with the remainder being removed by the apogee 
burn. In this conventional method, the final conditions at the end of each burn are known 
and the mass can easily be maximized progressively phase by phase if the second and 
third impulse sizes are unspecified. 

Now suppose that the total impulse of the third burn is fixed. Then the transfer 
orbit must be constructed such that the AV available from the third impulse is exactly 
that required to place the vehicle in circular synchronous equatorial orbit. If the AV 
available from a fixed third total impulse is less than that required to circularize and 
equatorialize the orbit for the mass available from a conventional ascent and second im- 
pulse, then the trajectory to transfer orbit insertion must be altered to reduce the AV 
required of the third impulse. This can be done by reducing the required plane change 
and the AV required for circularization. For reasons described at length in the RE- 
SULTS AND DISCUSSION section, the unconventional trajectory needed to reduce the AV 
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required of the apogee motor varies in many respects from the conventional profile. The 
most dramatic changes are a non circular parking orbit, nontrivial inclination reduction 
by the second impulse, and a significantly nonequatorial latitude for the second impulse. 
As is desired, the changes result in lowering the AV required of the fixed apogee motor. 
However, in this unconventional profile, the final conditions required at the end of the 
ascent and second impulse are unknown. They might be determined by varying those 
final conditions parametrically until the optimum is obtained. However, because of the 
number of variables, this process is clumsy and time consuming. 


Calculus of Variations Solution 

A calculus of variations formulation was used to maximize the payload to circular 
synchronous equatorial orbit without resorting to a parametric search. The optimiza- 
tion of the atmospheric portion of the trajectory is omitted from the variational analysis 
since the steering is severely constrained by aerodynamic loading and heating limitations. 
The analysis considers the problem from the point in the trajectory where the atmo- 
sphere can be neglected to insertion into the final orbit. In addition to optimizing the 
steering, the durations of any unspecified burns and coasts are optimized while main- 
taining the specified perigee radius of the parking orbit. The analysis is presented in 
appendix B. It is derived in three-dimensional rectangular coordinates in a manner 
similar to reference 4. The equations for optimum burn and coast duration are obtained 
from an analysis similar to that used in reference 3. It is necessary to extend the 
analysis to include an intermediate boundary condition which specifies the perigee radius 
of the parking orbit at the end of the ascent. Additionally, the oblate Earth model must 
be added to the variational analysis. The effect of oblateness is not negligible in tra- 
jectories to circular synchronous equatorial orbit. Trajectories to that orbit are long, 
minimally around 6 hours. Oblateness is the major perturbing force during most of a 
trajectory. Because of the large change in inclination required to perform the mission, 
any perturbation in the inclination, thus increasing or decreasing the amount of plane 
change required of the propulsion systems, affects the final mass and should be con- 
sidered in the analysis. 

The trajectories are numerically integrated to incorporate a nonimpulsive vehicle 
model and to include the effects of oblateness and small thrusts over long periods of time 
which cannot be conveniently treated impulsively. 

The analysis presented in appendix B requires the solution of a two-point boundary 
value problem. The solution to the two-point boundary value problem for the circular 
synchronous orbit problem with a fixed apogee burn and parking orbit coast time re- 
quires satisfaction of a minimum of eight final conditions with an equal number of initial 
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conditions. The number and specific initial and final conditions are explained in appen- 
dix B. 


Procedure 

A simple Newton -Raphson iteration scheme was used to solve the two-point boundary 
value problem. This scheme was used successfully with as many as 12 iteration var- 
iables. For further explanation of the iteration scheme, see reference 5. 

The partial derivatives required for the iteration scheme were obtained by integra- 
ting the adjoint equations. These were obtained as in reference 4. Solutions were ini- 
tially obtained using a spherical Earth model for the adjoint equations, but it was found 
that including the oblateness terms improved the convergence properties of the problems. 
In some problems of this type, it was found that including the oblateness terms was 
necessary to obtain convergence. 

It was difficult to obtain solutions to these problems because of the high degree of 
nonlinearity of many of the derivatives as well as the difficulty of guessing the initial 
values of the thrust angle in pitch and yaw and their rates. A technique was devised to 
systematically proceed from a simple, easily converged problem to the final solution. 
This technique is described at length in appendix C. Other techniques, such as gradient 
methods, might avoid some of the difficulties associated with the Newton -Raphson tech- 
nique. However, the method described in the appendix is convenient, straightforward, 
and adequate. After getting one solution, proceeding to others in the region of interest 
is not difficult. 


RESULTS AND DISCUSSION 

The Applications Technology Satellite (ATS)-E mission is a circular synchronous 
equatorial orbit mission. The ATS program has the objective of advancing technology 
in areas which may have application to future spacecraft. The experiments which are 
conducted are spacecraft, technology, and science oriented. 

The spacecraft-oriented experiments on the ATS-E provide information on power 
supply and control systems, a gravity -gradient stabilization system, resistojet and ion 
micropound thrusters, and synchronous environment. The technology-oriented experi- 
ments are primarily concerned with communication problems. The scientific-oriented 
experiments gather data on the particle (electron and proton) distribution and flux and 
the character of the electric and magnetic fields at synchronous altitude. 

The launch vehicle for the ATS-E mission was an Atlas -Centaur and the solid apogee 
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motor was a part of the spacecraft system. The apogee motor total impulse was sized 
for the early ATS missions on the Atlas -Agena launch vehicle, which has less payload 
capability than the Atlas -Centaur. A desire to increase the spacecraft mass led to the 
change to the Atlas -Centaur. The apogee motor, although smaller than optimum for the 
larger vehicle, remained unchanged. 

The Atlas -Centaur is a two-and-a-half stage vehicle. The Atlas is propelled by two 
booster engines and one sustainer engine. The booster engines are jettisoned at a pre- 
determined acceleration level. The sustainer engine continues to burn (sustainer solo). 
The Centaur insulation panels and then the payload fairing are jettisoned during this phase. 
The sustainer solo ends at propellant depletion and the Atlas stage is jettisoned. After 
about 10 seconds, the Centaur engines, burning hydrogen and oxygen, ignite and burn 
until the desired parking orbit is reached. During the parking orbit, a hydrogen peroxide 
propulsion system is used to maintain a very small acceleration for propellant retention 
and for attitude control. At the end of the parking orbit, the Centaur engines burn again 
until the proper transfer orbit is achieved. After engine shutdown, the Centaur control 
system acquires the proper orientation for the spacecraft burn, the Centaur and the 
spacecraft separate, and the spacecraft is spun up for stability. The spacecraft coasts 
up to the proper altitude maintaining the separation attitude. The spacecraft motor 
burns to place the spacecraft in the final orbit. The spacecraft apogee motor has a 
thrust of 22 240 newtons and an effective specific impulse of 279. 1 seconds. The total 
impulse available from the motor is 950 900 newton-seconds, which corresponds to a 
propellant load of 347 kilograms. 

The trajectory starts with a short vertical rise followed by a rapid pitchover phase 
in the desired azimuth direction. The amount of pitchover determines the amount of 
lofting during the atmospheric portion of the trajectory. The remainder of the atmo- 
spheric phase (which is assumed to end at booster stage jettison) is flown with a near- 
zero angle of attack steering program (described in ref. 5) to minimize vehicle heating 
and aerodynamic loads. The thrust direction is constrained to be parallel to the launch 
azimuth plane, which is established at launch. 

Since the ATS-E spacecraft motor has a fixed propellant load, the trajectory must 
be designed such that the AV required at apogee of the transfer orbit is exactly that 
required to place the spacecraft in the desired final orbit. As mentioned earlier, the 
ATS-E motor is much smaller than optimum. The Atlas-Centaur can put more mass in 
a conventional transfer orbit than the apogee motor can place in circular synchronous 
equatorial orbit. Therefore, an un conventional trajectory is required to lower the AV 
required of the apogee burn. 

An optimum unconventional trajectory was obtained for the ATS-E mission to cir- 
cular synchronous equatorial orbit. The AV required of the apogee motor is reduced 
by decreasing each of the two components which together make up the total AV - that 
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needed to circularize the orbit and to reduce the inclination to zero. The AV for cir- 
cularization is reduced by increasing the horizontal velocity at apogee of the transfer 
orbit without adding radial velocity. Any radial velocity would have to be removed by 
the apogee burn. Increasing the horizontal velocity at a fixed apogee radius is equiva- 
lent to raising the perigee radius of the transfer orbit - thereby decreasing the elliptic- 
ity of the transfer orbit. 

The AV required at apogee for reducing the inclination to zero is decreased by 
lowering the inclination of the transfer orbit. However, raising the velocity at apogee 
increases the AV required for inclination removal at a fixed transfer orbit inclination. 
Therefore, the combination of the two methods represents a compromise which is op- 
timized as part of the total problem. Most of the inclination reduction is performed by 
the second burn near the equator and only a small part is accomplished in the ascent to 
parking orbit. 

The characteristics of the optimum parking orbit are changed from the conventional 
profile to increase the perigee radius of the transfer orbit. An elliptical rather than 
circular parking orbit is used to raise the altitude of the second burn. The perigee 
radius of the optimum parking orbit remains limited by aerodynamic heating consider- 
ations at some acceptable value. Since injection into the parking orbit occurs near 
perigee, the vehicle must coast along the ellipse to a higher radius. Because of the 
limitation of the coast duration for the ATS-E mission, the second burn was required to 
occur near the first equator crossing. (From tracking or other considerations, a second 
(or greater) equator crossing could be chosen for the second burn, which would increase 
the parking orbit coast time by a half period (or more). ) The latitude of the second 
burn is no longer equatorial as in the conventional case since the optimum position for 
raising the perigee radius and decreasing the inclination is dependent on radius and 
velocity as well as latitude. The parking orbit coast time is greater for this unconven- 
tional profile since the time to the equator is greater for an elliptical than for a circular 
parking orbit and, additionally, the second burn occurs significantly south of the equa- 
tor. Optimum true anomalies are found for the beginning and end of the parking and 
transfer orbit coasts. In addition, the optimum combination of the changes just described 
as characterizing the unconventional profile is selected. 

The desired final inclination for the ATS-E mission is not exactly zero. The per- 
turbations of the Sun, Moon, and oblateness of the Earth cause a spacecraft to drift from 
an exactly equatorial orbit. Since zero inclination is not a stable condition, a final orbit 
inclination yielding the smallest average inclination over the lifetime of the spacecraft 
is desired. Small final inclinations with the proper inertial ascending node are found to 
yield acceptable inclination over the lifetime of the satellite. The particular combin- 
ations of final orbit inclination and ascending node are functions of the positions of the 
Sun and Moon, which are in turn functions of launch time and date. Therefore, data 
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were obtained for payload to circular synchronous orbits as a function of final inclin- 
ation. Negative inclinations are included in the data. This convention indicates that the 
node has been switched approximately 180° by the apogee burn. 

The Atlas-Centaur has a 25-minute limitation on parking orbit coast time for the 
mission. Therefore, inclusion of that constraint is necessary for realistic determi- 
nation of vehicle capability. However, optimizing the coast time provides a more dra- 
matic and obvious demonstration of the optimization procedure. 

Figure 3 presents separated spacecraft mass as a function of launch azimuth for 
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Figure 3. - Separated spacecraft mass as function of 
launch azimuth. 
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final inclinations of (-)2° and 5.25° for both optimum and 25-minute parking orbit 
coast times. Separated spacecraft mass is the mass of the spacecraft when it is sep- 
arated from the Centaur vehicle. This figure shows that the separated spacecraft mass 
is rather insensitive to launch azimuth. Hence, for simplicity, launch azimuth is fixed 
at 90° for the remaining figures. 

Figure 4 shows the separated spacecraft mass as a function of final inclination. The 
separated spacecraft mass decreases as final inclination decreases. Figures 5 and 6 
show the effect of final inclination on the transfer orbit inclination and inertial velocity at 
apogee. As might be expected, as the final inclination decreases, so does the transfer 
orbit inclination because of the nearly fixed amount of inclination change supplied by the 
small apogee motor. 

As might not be expected, the velocity at apogee also decreases as final inclination de- 
creases, especially for the case of optimum coast times. Figures 7 to 12 show why this 
occurs. Figure 7 shows the latitude of the second Centaur engine start as a function of 
final inclination. Since the second burn is required to remove more inclination as final 
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Figure 5. -Transfer orbit inclination as function of 
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Figure 6. - Inertial velocity at apogee motor ignition as 
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inclination decreases, it is advantageous to move the burn nearer the equator for more 
efficient plane change. Figure 8 shows that the longitude of second burn start also de- 
creases as final inclination decreases. These trends decrease the parking orbit coast 
arc as final inclination decreases. This is reflected in a decrease in the true anomaly 
at second Centaur cutoff, as seen in figure 9. Figures 10 to 12 show additional effects 
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Figure 9. - True anomaly at second main engine cutoff 
as function of final inclination. 
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Figure 10. - Parking orbit coast time as function of final 
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of moving the second burn nearer the equator. It decreases the parking orbit coast 
time, the altitude of the second burn, and the apogee altitude of the parking orbit. These 
all occur as a result of the decrease in parking orbit coast arc. These figures show why 
the apogee velocity is decreasing as final inclination decreases. The perigee radius of 
the transfer orbit decreases as the altitude of the second burn decreases. The apogee 
altitude of the transfer orbit is almost constant at synchronous altitude; hence, as peri- 
gee decreases, so does apogee velocity. 

Figures 5 and 6 also indicate that more AV is required of the apogee motor as 
final inclination decreases. It can be seen that both the plane change and circularization 
AV are increasing. However, figure 4 shows that the spacecraft ignition mass is de- 
creasing, which increases the AV capability of the apogee motor. 

Figure 13 shows the percentage of the Centaur propellant used in the first burn. The 
figure shows that as the final inclination increases, the first burn duration increases as 
the apogee altitude increases (fig. 12). 

The final longitude as a function of final inclination is shown in figure 14. The sat- 
ellite remains at the longitude indicated only when the inclination is zero, the orbit cir- 
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Figure 11. - Altitude at second main engine start as 
function of final inclination. 
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Figure 13. - Percentage of Centaur propellant burned 
during first burn as function of final inclination. 
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cular, and the altitude synchronous. For other inclinations, the position (latitude and 
longitude) of the satellite subpoint describes a figure eight on the surface of the rotating 
Earth. The longitudes indicated in figure 14 are injection longitudes, not necessarily 
the longitude at which the equator crossing occurs. For small inclinations, the longitude 
does not vary greatly during the period of the orbit. Figure 14 shows that longitude de- 
creases as final inclination increases. 

Now consider the limitation of parking orbit coast time. Twenty -five minutes is 
less than optimum for all the final inclinations considered, as seen in figure 10. The 
difference in separated spacecraft mass is shown in figure 4. As seen from these fig- 
ures, as the difference between the optimum and limited coast times decreases, the loss 
in payload due to coast time limitation decreases also. 

The coast time limitation reduces the advantage of raising the apogee of the parking 
orbit as final inclination increases. The energy required to raise apogee does not yield 
the payload increases available with optimum coast time since altitude can not be acquired 
as efficiently in the shorter coast time. The energy is better spent by the second burn 
to reduce the inclination of the transfer orbit. This is reflected in several of the fig- 
ures. In figure 5, the transfer orbit inclination for the coast limited case lies well be- 
low the optimum case. The lower second burn altitude is reflected in the lower velocity 
at apogee of the transfer orbit, as seen in figure 6. Because the parking orbit charac- 
teristics do not vary greatly with final inclination, the latitude and longitude of the second 
burn and the true anomaly at second Centaur cutoff are nearly constant. These may be 
seen in figures 7 to 9. The conclusions which may be drawn for the percentage of Cen- 
taur propellant used in the first burn and final longitude (figs. 13 and 14) are similar to 
those for the optimum coast case. 


Conventional Circular Synchronous Equatorial Orbit Trajectories 

The analysis and techniques developed to optimize the unconventional trajectory 
were applied to the conventional problem. References 1 and 2 show that the conventional 
trajectory profile for reaching a circular synchronous equatorial orbit is optimum if the 
second and third burns are impulsive and free for optimization. As mentioned in the 
INTRODUCTION, the referenced analyses assume that a circular orbit about a spherical 
Earth has already been achieved. The effects on the optimum profile of an oblate Earth, 
a nonimpulsive vehicle model, and the ascent to parking orbit were not considered. 

These effects are investigated with the present, more complex, analysis. 

The vehicle model chosen for the second study is the Atlas -Centaur, except that the 
Centaur stage of the vehicle is simplified so that the trajectory characteristics will not 
be obscured by vehicle peculiarities. In particular, the Centaur stage is used with 
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constant thrust and specific impulse from Atlas sustainer jettison through injection into 
circular synchronous equatorial orbit. No solid apogee motor is used. All jettisoned 
weights were eliminated in that period and only the main engines were used. The park- 
ing and transfer orbit coasts were without thrust and mass flow. 

The durations of the first, second, and third Centaur burns are optimized, as are 
the two coast phases. The resulting optimum trajectory had a circular parking orbit and 
the optimum launch azimuth was 91°. About 0. 3° of inclination was removed during the 
ascent to the circular parking orbit. Although a 91° launch azimuth is only slightly bet- 
ter than using a 90° launch azimuth, moving the trajectory closer to the equator more 
quickly makes plane change in the ascent more advantageous. The second Centaur bum 
spans the equator and lowers the inclination by 2. 1°. The remaining inclination is re- 
moved at apogee of the transfer orbit as the trajectory is circularized. 

The analysis used for this conventional problem is identical to that used in the un- 
conventional problem except for some modifications necessary because the optimum park- 
ing orbit is circular. The perigee altitude constraint becomes more complex for a cir- 
cular orbit. This is treated in appendix B. 

The results confirm that the conventional trajectory profile is essentially optimum 
even with an oblate Earth model and a nonimpulsive vehicle simulation. 


SUMMARY OF CONCLUSIONS 

Analysis and results are presented for trajectories to circular synchronous equa- 
torial orbit where the apogee motor is fixed at a smaller than optimum total impulse and 
for a conventional optimum configuration. The results for the small apogee motor case 
were obtained for the ATS-E mission, which used the Atlas-Centaur launch vehicle. 

The results confirm that the conventional trajectory profile is essentially optimum 
vehicle-apogee motor combinations where the apogee motor is smaller than optimum. 
More important, the results demonstrate that optimum trajectories to circular synchro- 
nous equatorial orbits may be obtained with detailed and hence complicated vehicle 
models for unconventional (small apogee motor) as well as conventional (optimum burn 
time) trajectory profiles. These results may be obtained without resorting to exotic 
mathematical procedures for solving the two -point boundary value problem. These re- 
sults were obtained with a simple Newton -Raphson iteration scheme. The partial der- 
ivatives were obtained by integrating the adjoint equations. The simple iteration scheme 
with the integrated partial derivatives is able to obtain solutions to the highly nonlinear 
two -point boundary value problem even when the number of initial and final conditions 
reaches twelve. 

The conventional circular synchronous equatorial orbit problem is also investigated 
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with the complex analysis used for the unconventional problem. The results demonstrate 
that the conventional trajectory is essentially optimum even with an oblate Earth model 
and a nonimpulsive vehicle simulation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1970, 

180-06. 
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APPENDIX A 
SYMBOLS 


c 

first integral of Euler- 

z 

unit vector pointing at north 


Lagrange equations, kg/sec 


pole 

E 

energy per unit mass, m^/sec^ 

|3 

mass flow rate, kg/sec 

e 

eccentricity, dimensionless 

e 

jump factor 

F 

functional defined by eq. (B4) , 

V 

Lagrange multiplier, kg/sec 


kg/sec 

X 

Lagrange multiplier, kg-sec/i 

f 

unit thrust direction 

M 

Lagrange multiplier, kg/m 

G 

gravity acceleration, m/sec^ 

a 

Lagrange multiplier, dimen- 

G m 

spherical Earth gravity con- 

3 2 

stant, mysec 

(P 

sionless 

yaw attitude, deg 

G 1’ G 2 

components of oblate gravity 
acceleration, m/sec^ 

* 

pitch attitude, deg 


Superscripts: 

g 

intermediate boundary equa- 



tion 

f 

final 

h 

angular momentum per unit 

o 

initial 


mass, rn^/sec 

- 

time derivative 

J 

functional to be minimized, 

- 

vector 


kg 

A. 

unit vector 

m 

mass, kg 

Subscripts: 

N 

total number of stages 

d 

desired 

P 

semilatus rectum, m 

f 

final 

r 

radius, m 

i> \ 

Z,m,n J 

stage numbers 

r p 

perigee radius, m 


S 

variational switching function, 

0 

initial 


dimensionless 

pk 

parking orbit 

T 

thrust, N 



t 

time, sec 



V 

velocity, m/sec 
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APPENDIX B 


DERIVATION OF OPTIMUM CONTROL 

As mentioned in the ANALYSIS section, the optimization of a trajectory to a cir- 
cular synchronous equatorial orbit may be considered as the problem of optimizing a 
multistage launch vehicle to a particular final orbit. The optimization problem to be 
considered here begins at booster jettison, which is assumed to be at a fixed position and 
velocity. The sustainer portion of the Atlas continues until propellant depletion. The 
sustainer is jettisoned and a few seconds later the first Centaur burn begins. Its dura- 
tion is variable and must be optimized. The perigee radius of the parking orbit which 
follows is fixed. The duration of the parking orbit may or may not be optimized. The 
parking orbit is not a true coast since a small acceleration is maintained for propellant 
retention. The duration of the second Centaur burn must be optimized, followed by an 
optimum transfer orbit coast (a true coast) and a final burn of fixed total impulse. The 
analysis presented in this appendix to solve this problem is a special case of the analysis 
in reference 3, with the additional constraint of a parking orbit perigee radius. 

The variational problem to be solved is to find the steering program and various 
stage durations which maximize the payload capability of a multistage launch vehicle to 
a specified final orbit. The trajectory must satisfy certain initial, final, and interme- 
diate conditions on the state variables. The thrust, propellant flow rate, and jettison 
weight for each stage are assumed to be constant. The equations of motion and con- 
straints for each stage may be written as 


_ T. 

v - G(r) - _i f = 0 

(Bla) 

m 


!o 

II 

l> 

(Bib) 

m + = 0 

(Blc) 

a a 

f • f - 1 = 0 

(Bid) 


where f is the unit thrust direction and G(r) is the oblate Earth gravity acceleration 
(ref. 6) , which may also be written 

G(r) = G^r, r • z)r + G 2 (r, r ■ z)z (B2) 
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Suppose that each stage of the vehicle is numbered consecutively starting with the booster. 
For analysis purposes a stage change occurs when the thrust and/or propellant flow rate 
changes and/or a mass is jettisoned. A Bolza formulation of the variational problem is 
used (ref. 7), and the functional to be minimized is written as in reference 3 as 

J = -m 



where the functional F. for each stage is 


Fi = X*[v-G- — f)+M‘ (r - v) + a(m + /3.) + rj(i • f - 1) 
1 1 m 1 1 


The resulting Euler -Lagrange equations are 


X + M = 0 


. G. X 

M + — — + (X • rJV-Gj - G 1 (X • r) - + (X • z)V-G 2 = 0 


• T i - - 

d L \ ■ f = o 


m 


(B5a) 

(B5b) 


(B5c) 


T. _ 

2 v f - -1 X = 0 (B5d) 

m 


A 

The optimum thrust direction f is obtained by combining equations (Bid) and (B5d) and 
using the Weierstrass E-test. This procedure results in 

f = X (B6) 

Integrals of the Motion 

Since F does not explicitly depend on time, an integral of the motion is 
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(B7) 




_ _ _ _ T. 

(L + X • G + p • v + — X - a/3. = 0 
1 m 

_ _ Q 

When a spherical gravity model is assumed (i.e. , G(r) = (G^/r)r), three additional 
integrals of the motion exist which are given by 


Xxv+pxr= constant 

Since X, p, r, and v are all continuous except where an intermediate boundary con- 
dition is imposed (as will be shown later), the three integrals are constant across staging 
points where continuity holds. However, for the oblate gravity model used in this 
analysis, only a single component of the previous vector integral is constant, as can be 
verified by differentiation with respect to time: 

(X X v + M x r) • z = constant (B8) 


Transversality Equation 


The transversality equation for this problem is 


dj 


N _ t. 

= y (C dt + X • dv + P • dr + a dm)^ 1 - dm^ 

i=2 i_1 


(B9) 


which is set equal to zero for an optimal solution. Reference 3 shows that X and p are 
continuous everywhere if there are no intermediate boundary conditions. If the inter- 
mediate boundary condition (assumed to occur at a staging point) is expressed as 


g(r,v) = 0 


(BIO) 


reference 8 shows that the discontinuities in X and p" are eV-g and eV^g, respec- 
tively. The variable e is used as an initial condition in the two -point boundary value 
problem to satisfy the intermediate boundary condition (eq. (BIO)). 

The equations that must be satisfied to optimize the duration of the powered and 
coast stages are derived in reference 3. The applicable results are presented here. Let 
j be the first optimized powered stage. Then for constant jettison weight the equation 
for optimizing stage l is 
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(Bll) 



= 0 


where o and f refer to initial and final values and the S functions are defined as 



a = 


- ■ G + M • v + — - x| 


*i 


for /3j ^ 0 


(B12a) 


S i s 0 for /3. = 0 


(Bl2b) 


The right side of equation (B12a) is obtained by using equation (B7). For coasting stages 
(0. = T. = 0) to be optimized, the equation 

C. = -(X-G+ju-v) = 0 (B13) 


must be satisfied for maximum payload. 

For free initial or final state variable x, the required or final condition for maxi- 
mum payload (ref. 4) is 


A • — + p. • — = 0 (B14) 

0x 9x 

Initial Conditions 

If the initial position and velocity are specified, the initial values of any five of the 
six X and may be used as variable initial conditions to satisfy the required final 
conditions of the two-point boundary value problem. To eliminate the difficulty associ- 
ated with guessing at values of the multipliers, the values of A and JT can be expressed 
in terms of pitch and yaw attitude (\p and cp) and rates and cp). These equations 
may be found in appendix C of reference 4. The values of A and /I are then calculated 
from 


X - XX 


(B15a) 


M = - AA — AA 


(B 1 5b) 
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The value of X can be set equal to unity without loss of generality. The initial value of 
X can be calculated in closed form, as will be shown by the following development. 


Final Conditions 

Final conditions for both the conventional and unconventional synchronous equatorial 
orbit mission require a circular orbit at synchronous orbit altitude with prescribed in- 
clination. If the required inclination is nonzero, both the longitude of the ascending node 
and the injection point in the final orbit are free for optimization. As shown in refer- 
ence 4, the corresponding auxiliary variational final conditions are 

(xxv"+flxr) • z = 0 (B16a) 


and 


(XX V + 11 x r) • (r x v) = 0 (B16b) 

If the desired inclination is zero, equations (B16a) and (B16b) degenerate into one equa- 
tion (zero inclination is equivalent to two final conditions, r • z = 0 and v • z = 0), and 
only equation (B16a) must be satisfied. 

Since equation (B16a) is a constant of the motion (eq. (B8)), it may be satisfied at 
the beginning of the trajectory and used to calculate X. However, it must first be veri- 
fied that jump discontinuities in X and "p at intermediate boundary points do not change 
the value of the constant. This requires that 

(Vyg X v + V-g X r) • z — 0 (B17) 


It will be shown later that equation (B17) is satisfied for all functions g used herein. 

The calculation of X proceeds as follows: 

(Xxv + juxr)' z = (XX x v) • z - (XX x r) • z - (XX x r) • z = 0 

X = MX X V - X X r) • z (B18) 

(X x r) • z 

Computing X with equation (B18) guarantees that equation (B16a) will be satisfied. 


21 


I 



Intermediate Conditions 


As explained earlier, it is necessary to constrain the perigee radius at injection into 
the first parking orbit. Otherwise, the optimum solution would result in the parking 
orbit injection and/or the equator crossing occurring at very low altitudes, thus violating 
spacecraft heating constraints. Therefore, the intermediate constraint is 

g(r,v) = r p - r p d = 0 (B19) 

where the desired value corresponds to the perigee altitude. If equations found in ref- 
erence 9 are used, equation (B19) can be written as 

- r n . = 0 (B20) 

1 + e p ’ 


where 


i_ u 

p = n ' n (semilatus rectum) 
G m 


(B2 la) 


e = 


r 


(eccentricity) 


g; 


m 


(B21b) 


- - G 

E = v ‘ v — — (energy per unit mass) 

2 r 


h = r x v (angular momentum per unit mass) 


(B21c) 

(B21d) 


The required gradients are calculated to be 




h(h x r) - r^v 


eG; 


m 


(B22a) 


(vxh) - r -S 


V F g = 


m 


(B22b) 
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It is easily shown that equation (B17) is satisfied for the gradients in equations (B22a) 
and (B22b). In fact, equation (B17) is satisfied for any function g of r, v, h, and 
r • v. For such a function g, 


V-g = -Mv + — (hXr)+ — r 

3v 3h 3(r • v) 

V_g = + ^(vxh)+ — v 

r 9r 0h 3(r • v) 


and 


V— g x v + V^g Xr = iivXv + -^(hX r) x v + — ^ r x v 

v 3v 3h 3(r • v) 


+ ^£rXr + — (vxh)x r + — vXr=0 

3r 3h 3(r • v) 


Hence the value of X x v + p. x r is unaffected by the jump in X and p. 


Conventional Case 

As discussed in the RESULTS AND DISCUSSION section the circular parking orbit 
problem required additional analysis. For the conventional synchronous orbit case, it 
became obvious during the iteration that the trajectory was converging to a circular 
parking orbit. However, as the parking orbit became more and more nearly circular, 
it became increasingly more difficult to obtain convergence. The problem was traced 
to the gradients of perigee radius (eqs. (B22)). It is seen that the gradient becomes in- 
determinate as the eccentricity approaches zero. Attempts at resolving this indeter- 
minacy showed that the indeterminacy is fundamental; that is, the gradient direction and 
magnitude are not defined for a circular orbit. 

In order to resolve the indeterminacy and obtain convergence, the constraint was 
restated as requiring a circular orbit at the desired radius. This requires three con- 
straints: 


gl = r - r d = ° 

(B23a) 

g 2 = v - v d = ° 

(B23b) 
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(B23c) 


g 3 = r • v = 0 

The corresponding gradients are 

V^gj = r, V-g x = 0 (B24a) 

V^g 2 = <5, V v g 2 = v (B24b) 

V-g 3 = v, V-g 3 = r (B24c) 

Three jump scale factors are required for this case in order to converge on the 
intermediate conditions (eqs. (B24)). With this method, convergence was easily ob- 
tained. 

It was verified that the optimum conventional case requires a circular parking orbit. 
This was accomplished by comparing the payload for a circular parking orbit case with 
the payloads for slightly noncircular parking orbits whose perigees are equal to the cir- 
cular parking orbit altitude. A map of payload points was obtained which indicates that 
a circular parking orbit is optimum for this problem. 


Boundary Value Problem 

For the ATS-E mission, both fixed and optimum parking orbit coast times were con- 
sidered. The transfer orbit coast time was always optimized, however, along with the 
durations of the first and second Centaur burns. Based on the preceding discussion of 
the transversality equation, the initial and final conditions for the two-point boundary 
value problem are shown in table I for the case where the parking orbit coast time was 
optimized. 

If the desired final inclination is nonzero, then (r • z) and (v • z) = 0 in table I are 

replaced by i H and (JI x r + Xxv)-(rxv) = 0. If the parking orbit coast time is fixed, 

k- 1 

then an initial and final condition are removed. These are t^ and ^ ] (s? - S? + jj = 0. 

i=j ' ' 

It should be recognized that there may be any number of fixed stages between t. and t. , 

J ^ 

etc. Also, the last three final conditions are evaluated at intermediate points in the 
trajectory. 

For the ideal synchronous orbit mission, the two coast phases, as well as the three 
Centaur burns, were all optimized. The resulting initial and final conditions in the two- 
point boundary value problems are given in table n. 
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TABLE n. - INITIAL AND FINAL CONDITIONS 


TABLE I. - INITIAL AND FINAL CONDITIONS 
FOR UNCONVENTIONAL CASE 


Initial conditions 

Final conditions 

p 

E d 

p 

r a 

<P 


<P 

e 

r . (parking orbit) 
p , a 

(r * z) = 0 

t. (first Centaur burn) 

o 

ii 

<N 

t^ (parking orbit coast) 

E (4 - s r + i) = 0 

H 

tj (second Centaur burn) 

E (4 - s ° + i) = 0 

i=3 

t m (transfer orbit coast) 

(A • Tj + TZ • v) = 0 


FOR CONVENTIONAL. CASE 


Initial conditions 

Final conditions 

P 

E d 

P 

r d 

<P 

r d 

V 

r pk, d 

e l 

v pk, d 

e 2 

(?■ V) pk = 0 

€ 3 

(r • z) = 0 

t. (first Centaur burn) 

(v • z) = 0 

t k (parking orbit coast) 

(F • v + X • G) k = 0 


l-l / 

t^ (second Centaur burn) 

E (4 - s h ) = 0 


i=j ' 

t m (transfer orbit coast) 

(M • v + A • G) m = 0 


n-1 / 

t (third Centaur burn) 

E (4 - 4 + i) = ° 


i=j ' 


As in the ATS-E case, if the final orbit inclination is nonzero then (r • z) = 0 and 
(v • z) = 0 in table n are replaced by i^ and (]Tx r + Xx v) • (r x v) = 0 . 
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APPENDIX C 


TWO-POINT BOUNDARY VALUE PROBLEM 

The following technique was devised to systematically proceed from a simple, easily 
converged problem to the solution of the two-point boundary value problem for a circular 
synchronous equatorial orbit. 

A trajectory is obtained to a slightly elliptical (parking) orbit with the desired peri- 
gee radius without plane change with a 90° launch azimuth. This problem converges 
easily. Then the ascent burn time is fixed at the value obtained and a variable length 
parking orbit coast, a fixed parking orbit perigee radius and second burn are added. 

This problem is targeted to the desired apogee and 180° argument of perigee for first 
equator crossing second burn. An inclination decrease of about 2° is then added to these 
final conditions and the problem is retargeted to the augmented final conditions. Now 
the transfer orbit coast (variable) and the apogee burn (fixed or variable) are added. 

This trajectory is integrated to the end with the converged initial guesses from the last 
step. The final conditions achieved will frequently be far from a circular synchronous 
equatorial orbit. However, specify the final conditions actually achieved as the desired 
ones, and optimize the problem. The parking orbit coast, second burn, and transfer 
orbit coast durations will change. Now alter the achieved final conditions toward the 
desired ones judiciously in steps, retargeting at each step. In this manner, the desired 
final orbit conditions may be obtained. Now the ascent bum duration may be optimized. 

If during the process of optimizing the ascent the parking orbit moves toward being cir- 
cular, a fundamental difficulty is encountered in the equations for optimizing the problem 
with the perigee constraint. This is discussed further in appendix B and the RESULTS 
AND DISCUSSION section. Any sizable change in a constraint or final condition is best 
achieved by proceeding in steps. The problem is quite nonlinear. Attempts to obtain 
initial conditions as functions of the final conditions by extrapolation of data obtained 
from converged problems were made. They were generally unsuccessful. 
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